Next-nearest-neighbor spin-spin and chiral-spin correlation functions in generalized 

XXX chain 
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We develop a simple technique for calculation of next to nearest neighbor spin-spin and chiral- 
spin correlation functions in inhomogeneous XXX model. Exact expression of the chiral-spin order 
parameter as a function of the model parameter, lj, is analytically found. Using the same method 
we also calculate the next to nearest neighbor spin-spin correlation function. In the limit a; — > it 
reproduces the known result for the vacuum expectation value of the next to nearest neighbor spins 
in the standard Heisenberg spin chain. The technique is simple and can be extended for calculation 
of next to next to nearest neighbor correlation functions as well as for calculation of correlation 
functions in XXZ model. 

PACS numbers: 71.27.-|-a, 71.10.-w, 71.10.Hf 



I. INTRODUCTION 

Presently, calculation of correlation functions in 
strongly correlated electron systems is one of the most 
important tasks in physics of low dimensions. Corre- 
lation functions of spins at large distances are directly 
related to the observable quantities. They define a set 
of critical indices identifying the universality classes of 
different phases^*^'^'^. Notably, correlation functions of 
spins in the homogeneous Heisenberg spin chain (XXX 
model) at short distances are not directly related to the 
universality class of the phase, but, as it appeared^, they 
play an important role in the perturbative investigations 
of ladder models, that are not integrable. An important 
example is the Haldane ladder model^, that can be ap- 
proached as an Heisenberg spin chain perturbed with the 
chiral-spin order operator, i7„((t„+i x an+2), (which was 
considered earlier in connection with spin-liquid ordered 
phase in Ref0) and by the product of next to nearest 
neighbor spins,- ((t„i7„+2). In order to investigate the 
free energy of this long term attractive model and ana- 
lyze the phase space, one needs to know the correlation 
function of the next to nearest neighbor spins and the 
chiral-spin order parameter Xn = (f?n(CTn+i x (Tn+2))- 

The aim of the present work is the calculation of 
abovementioned correlation functions for a generalized 
XXX model^'-ii^iiiii^iii^ by a new simple technique. This 
model is defined and studied in Section II. The stud- 
ies of Section II involve analysis of the Hamiltonian(s) 
and other conserved currents, construction of the Trans- 
fer matrix and obtaining the set of Bethe equations that 
describe the energy spectrum. Two particular cases of 
the generalized XXX model are the ladder systems with 
N sites at each chain given by Hamiltonian operators 
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2JV 



n=l 



H = 2^2 [^n^n+l - 1] + ^ [an(Jn+2 - 1] 

n=l 

2N 

+t^^(-l)"o'„_l((T„ X CT„+l), (1) 



and 



2N 



?^' = c^2^(-i)"[a„a,r+2-i] 

2N 

-Uj'^an-liSn X 0^„+i). (2) 



As a manifestation of the effective workability of the de- 
veloped technique in Section HI we first calculate the 
exact correlation functions of the next to nearest neigh- 
bor spins, ^(a;) = {unOn+2 — I), as a function of to, which 
for the Hamiltonian Eq. ([1]) has the following asymptotic 
behavior: 

-161og(2)-H9C(3)-I,41u;2, 
. for cj < 1 

e(-) = <[ _2,77+i^, 

for cj > 1. 

(3) 

Here C,(x) is the Riemann zeta function. In the limit 
= 0, we reproduce the known result for the expectation 
value of next to nearest neighbor exchange operator in 
the standard XXX chain, first obtained by involving the 
Hubbard model^**. 

More recently, within the general approach of multiple 
integral representation of correlation functions, that was 
formulated in Ref. ,15i and investigated further in Refs.flGl 
and [13, correlation functions {undn+z) and {unOn+i) 
were evaluatecU^ii^. These functions were evaluated for 
XXX model in zero magnetic field. One of the advan- 
tages of our technique is the possibility to include also 
external magnetic field. This can be done by involving 
the magnetic field into the integral equations for the den- 
sities. In some cases such integral equations can be solved 
by approximate methods^. 

Our main result, however, is the analytical calculation 
of the chiral-spin order parameter, ((7„((t„+i x an+2)) = 
(— I)"x(^)j as a function of staggering parameter, w. The 
full expression and the derivation of x„ are presented in 
Section III. Here we present only the asymptotes of the 
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function x('^)j 

r [87 + 8V'(l/2)-V"(l/2)+^"(l)]c., 
x(w)= < for w <C 1, 

i^, foru;»l, 

(4) 

where 7 is the Euler constant, the function ip is the 
Digamma function and the coefficient of hnearity for 
small w is « 3, 33. 

We would like to note that this correlation function 
was previously investigated numerically^^ and the match 
between curves plotted from our exact result (Fig. 2) and 
from the numerical simulations is perfect. 



The interesting feature of the Transfer matrix 0(A; w), 
given by Eq. ([5]), is that instead of usual XXX Hamil- 
tonian this yields two different local Hamiltonian opera- 
tors, Hi and H2, that are proportional to the logarith- 
mic derivative of t at two different points: X = i/2 and 
A = uj + i/2. Respectively, their explicit forms are 

Hi = 2i{l + uj'^)dxhiT\^^.^^- N{2 + uj^ -2iuj) 

2N N 

= ^ [CTnCTn+l - 1] + X! [^■2k<y2k+2 - 1] 

n=l k=l 

— ^f?2fe((?2fc+l X '?2fe+2)} , (8) 



II. INHOMOGENEOUS CHAIN AND A FAMILY 
OF COMMUTING OPERATORS 

In one dimension, there is a variety of quantum exactly 
solvable models of interacting spins. Some of these mod- 
els involve interactions between nearest neighbor spins 
and also spins that are far from each other. Usually 
the further neighbor interactions and other non-localities 
come from additional anisotropy parameters. We will 
consider a family of models with nearest neighbor, next 
to nearest neighbor and triangular (zig-zag) interactions, 
which stem from the Transfer matrix with the shift of 
the spectral parameters at each second sites: 

e(A;u;) = ^2Ar.Q(A)L2Ar_i^Q(A - w) • ■ • (5) 
■ ■ •^2,a(A)Li^Q(A - LO). 

Here La,fc(A) obeys the Rational Yang-Baxter relations 

^ai,a2(A - Ai)i„,ai(A)i„,a2(M) (6) 

with rational 

in,a(A) = (A- V2)I„,a + iPn,a, (7) 
i?a,6(A) = AIa,b+iPa,b. 

The permutation operator is given in terms of Pauli ma- 
trices as 

a 

With this construction, one has a commuting one 
parametrical family of transfer matrices r(A;a;) = 
tr,e(A;a;) : 

[T(A;c^),T(M;a;)] =0. 

This is a well known picture, see e.g.,— >^ii^, while we want 
to apply this in a slightly different contents. In order to 
be shorter we will miss some proofs and refer the reader 
to the nice review by L. D. Faddeevi^. 



i/2 = 2i(l + w2)9AlnT|^^./2+^ -iV(2-Hw2_^2ia;) 

2N N 

= ^ [(TnO-n+i - 1] + X! {^^ [^2k-l(72k+l - 1] 

n=l fe=l 
-I- W(T2/e-l(CT2fe X (T2/C+1)} . (9) 

These operators are commuting as they belong to the 
same commuting family. It is straightforward that the 
Hamiltonian operators H = Hi — H2 and H' = Hi + H2 
are exactly diagonalizable in the same framework. Their 
explicit forms are given by Eqs. H]) and 

Let quasi- shift operators be the monodromy matrices 
at the points A = i/2 and A = i/2 + uj: 

U+ = traQ{i/2-oj), (10) 
U- = traQ{i/2 + Lu;Lu). 

They obey the relation 

il+uj^)^V^ ^U+U^ = e'P, (11) 

where 

V ^ Pl.2P2,3---P2N-l,2N (12) 

is a shift on one unit n — > n -|- 1, and P is the physical 
momentum, which governs the shift n — > n + 2. Being 
defined in this way, the quasi-shift operators commute 
with the whole family of transfer matrices t(/x; w). 

Derivation of the Bethe Ansatz Equations (BAE) for 
t{ij,;uj) is standard, see e.g., Ref4. Starting from the 
reference state with all 2N spins up one can generate the 
eigenvectors of the transfer matrix in the sector with M 
overturned spins, parametrized by M complex rapidities 
A„ which obey BAE 

nXn + i/2){Xn-iu + i/2) Y ^ ^ A„ -Xk + i 

\{Xn-t/2) {Xn -LO-i/2)) 4" " " ~ * ' 

(13) 
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The eigenvalues of the transfer-matrix t[^]lj) have the 
form 



M 



Z A 11 Aji 

n— 1 



M 



n 



^ - A„ + i 



n— 1 



(14) 



This gives the quasiparticle momentum in the form 



M 

n 



(A„ +^/2)(A„ + V2) 



(15) 



and eigenvahies of Hi and i?2 as follows: 
i.,({A},.)^-2(l+.^)E^, 



(16) 



£^2({A},CJ) 



-2(1+-^)E 



1 



(A„-cc>)2 + l/4- 



The corresponding eigenenergies of Eqs. H]) and ^ are 
i?i + i?2 and El — E2 respectively. Now the picture is in 
some sense complete and one is in position to infer the 
thermodynamics of these models, based on the BAE and 
the energy relations. The particular Hamiltonian Eq. ([T]) 
was introduced in Ref. [ill and analyzed in detail o^^'^^ . 
It has a singlet ground state with massless excitations. 
By involving a magnetic field with the Zeeman coupling, 
the system undergoes two phase transitions; two critical 
phases with different universality classes are discussed in 
Ref. lU The XXZ generalization of the model ^ was 
defined and investigated 



III. CHIRAL-SPIN AND OTHER 
CORRELATION FUNCTIONS 

It is well known that the expectation values of oper- 
ators not commuting with the Hamiltonian are not eas- 
ily accessible within the framework of Bethe Ansatz. In 
the case of inhomogeneous chain under consideration we 
have the additional parameter w, which breaks the trans- 
lational invariance n — > n -I- 1 and gives a possibility to 
calculate some simplest expectation values, that are valid 
also in the limit — > 0, corresponding to the well known 
XXX case. For our purposes we need the eigenvalues of 
quasi-translation operators, which follow from p4p and 
the definition PU)) : 



= (-1)^(1 + 



A„ + i/2 
A„ - V2' 



(17) 



A„ — uj + i/2 



X„ — Lu — i/2 

One can differentiate Eq. (fTOj) with respect to uj and 
let the relations 



N 



EiP2k,2k-l +^ jj 
T— ^ 



fe=l 



1-f Cj2 



N 



k=l 



2/c,2fc+l + ^ 



(18) 



N 

E 

fe=i 



iP2k,2k+l — 



N 



^ " ^ 1 -ha;2 



k=l 



(19) 



Let us use the relation, that is always valid when one 
has a parameter-dependent operator 0(r?) with the spec- 
trum On{r]) and normalized eigenstates \n) 



{n\d^d{r])\n) = 8,^0^(1]). 



(20) 



Then we find that 



{{^i}\anan+i\{fi}) 
2 

=F— i(l + u;'^)d^ ln[u±(w)] + 2iuj ~ 1 



M 



l-2ii±^^ 



N 



ni — 1 ' 



1/4' 



(21) 



where {/i} is any set of BAE solution. Given in the above 
form, it is readily calculable. It does not depend on n, 
even on the parity of n and shows that the model Eq. ([1]) 
can't have a dimerized phase. 

The next expectation value that we are going to eval- 
uate is the next to nearest neighbor (NNN) exchange, 
({/j,}|(T„(T„+2|{/J,}). For this purpose we divide the two 
Hamiltonian operators Eq. ([8]) by w and apply Eq. ([20)l . 
Subtracting contributions of nearest neighbor terms with 
the use of Eq. (I^T]) we obtain for even or odd site numbers 
the following relations: 
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(MI'?2fca2.+2 - 1|M) = ^d^^^^^^^^ + ^m\a^a^+i-im) (22) 

I\ LU UJ 



E 



and 



(Ml?2.-i5^2.+i - = ^ac. ^'^^^^'""^ + - (23) 

iV UJ 



M 



'^^A M^ + i 2C.2 (^„_^)2+l ^ ((^„-c.)2 + i)2 



It is yet unknown how to perform such summations in the 
case of finite N and M analytically. Instead we can eval- 
uate the sums in the important case of thermodynamic 
limit. 



N 



M 



M 

— = const, 



(24) 



when solutions of Eq. (|T3|) form bound states called 
strings'^ . In our rational case strings with arbitrary 
length n are possible. Consider the case where one has 
M„ bound states of ?i-strings 



1,..,M„ 



j = 1, 



(25) 



with real parameters of centers AJJ. Then one can take 
the product of n BAE for the same n-string and obtain 
an equation for real centers. The logarithm of these equa- 
tions gives 



0(A2/n) + 0((AS-c.)/n) 



— J 

N " 



iV 



(m,/3)7^:(n,a) 



(26) 



6nm(A^ - A™), 



K:' = AS + -(n + l-2j) + 0(cxp(-|<5|7V)), 



where B(\) = 2tan-i[2A] and 



0n7Ti(A) 



^(f^) + 2^(f^) + - + 2^(;i+^) + ^(;rT;;r) 



r 



for n ^ m, 
for n ~ m. 



(27) 



12 is an integer (half-odd integer) if N — Mn is odd (even) n-string /5„(A); the number of strings and holes between 



and satisfies 



t 



771 — 1 

~ 2Min(n^ m) — dnm- 



(28) 



A and A -1- dA is pn{X)NdX and pn{X)NdX, respectively. 
From Eqs. (|26p one obtains a system of integral equations 

a„(A) + Unix - Lu) ^ pn{X) + Pn{X) (29) 

+ E / ^nm(A ~ Ai)/9m(Ai)d/i, 



In the thermodynamic limit ()24|) , it is convenient to de- 
fine distribution functions of n-strings Pn(A) and holes of where Tnm{X) is a function defined by 

I 



flIn-'m|(A) -I- 2a|„_m|+2('^) + '^0.\n-m,\+i{^) + •■• 

(A) + a„+m(A) for n ^ m, 

, 2a2(A) + 2a4(A) + ... + 2a2„-2(A) -I- a2„(A) for n ^ m. 



(30) 
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and a„(A) is a function defined by 
a„(A) = -- 



In order to describe uj- derivatives of in the ther- 
modynamic hmit, we introduce a new function, Fn, as 



hm d^XlF,,iX,Lj). 



(31) 



this function, one can differentiate ((26)) with respect to 
UJ and use ([^5]) . In this way one finds: 

an{X-uj) - F„(A)[p„(A)+p„(A)] (32) 

/•OO 



For briefness we will miss the explicit w- dependence of Now we can rewrite the R.H.S. sums in (j^H) in terms of 
F. In order to find a characteristic integral equation for integrals: 

I 



and 



oo ^ 



n=l • 



—an{fl)Fn{li) 



2a;2 



-a„(^) 



1 + 

2tj 



n=l •' 



an{p)Fn{n) 



CJ^ — 1 1 + ti)^ , 

-a„(/i - w) + — a„(/i - uj){Fn{p.) - 1) 



2a;2 



2w 



(33) 



(34) 



Evaluation of the expectation value of triple interaction 
terms, the chiral-spin order parameter, can be done in a 
similar way. The answer is 



{{p}\^2k{$2k+l X l?2fe+2)|{M}) 
n—l ^ 



anip.)Fn{p) 



-aM + {1 + uj^)a'MF^{ti) 



and 



{{lJ'}W2k~l{^2k X ^2k+l\{fJ-}) 

1 



(35) 



(36) 



= 47r^ / pnip)dfiU- 

n=l •' 



an{p)F.a(p) 



an{p - c^) + (1 + LU^)a'Jii ~ u;)iFn{ti) - 1) 




FIG. 1: Next to nearest neighbor correlation function ^(o;) 
versus inhomogeneity parameter ui. 



Its solution 

Pi (A) 



1 



1 



2cosh7rA 2 cosh7r(A — w) 
and the corresponding solution to (j32p 

1 



Fi(A)pi(A) = 



For definiteness, let us evaluate the NNN expectation 
value for the ground state of the Hamiltonian Eq. ([T]). 
For this state, the densities are found to be zero for all 
the n- strings with n = 2,3, .. and l-holes^^. The system 
Eq. (123) reduces to the following simple integral equation for the NNN expectation value 
for n = 1: 



2 cosh7r(A — lo) 



(38) 



(39) 



can be found by the Fourier transform. The integrals in 
([55)1 can be easily transformed to the following expression 



ai(A)+ai(A-a;) = /9i(A)+ / Tii{X- p)pi{p)dp.. (37) 



f(w) = ({/i}|'?n'?n+2 - 1|{P}) 



(40) 



-lieu) 



1/(0) 



'-d^I{uj), 
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where /(w) is the integral which can be expressed via 
Digamma functions, Tp, as 



lico) 



dx 



-oo (^^ + 1/4) cosh7r(x — Lu) 



(41) 



We see that though the translational invariance n — > 
n+1 is broken, the vacuum expectation value of the NNN 
exchange operator does not depend on the site parity. 
Up to the sign factor, the same is valid for the triple 
interaction terms: (CT2fc('?2fe+i x (T2fc+2) = -{^2k-ii(^2k x 

In the limit a; — > 0, from Eq. (|40p . one will recover the 
known result for the expectation value of NNN exchange 
operator of Hcisenbcrg XXX isotropic chain, 



('7„(?«+2) = l-161og(2) + 9C(3). 



(42) 



This was calculated from the ground state energy of the 
Hubbard model in Ref . 14. We present the function ^ {us) 
in Fig.l. 

The chiral-spin order parameter ({^J■}\ffn{(^n+l x 
'?n+2)|{/^}} is also important, as it defines the measure of 
chirality of the state. Substituting the densities p and F 
for the ground state, Eqs. and into Eq. 
we obtain 



X{UJ) = (-l)"(a„(CT„+l X (T„+2)) 



(43) 



;/(0) 



2 + 4w^ 



I{Lu)-il+LU^)d^I{uj) 



This function is plotted in Fig. 2. We see a perfect 
match between our plot and the numerical simulations of 
Ref. [H. 

In conclusion, let us briefly comment on the possibil- 
ity of extension of the developed method to the third 
neighbor correlation functions, in particular of the type, 
((7„f7„+3). For this case one has to increase the level of 



inhomogeneity of the model by introducing two different 
shifts of the spectral parameter and consider the follow- 
ing monodromy matrix 

Q{X]UJi,UJ2) = L3Ar,a(A)i3Ar_i^a(A - t^i)£3A'-2,a(A - W2) 

••• i3,a(A)L2.a(A-cc;i)Li,,(A-cc;2), (44) 

which is defined on the lattice with 3A^ sites. With this 
construction, one again has an integrable model with 
commuting family of transfer matrices, but, contrary 
to case considered above, we will have now three local 
Hamiltonian operators. It is straightforward to derive 
their explicit forms, which are rather cumbersome and 
we don't bring them here. These operators contain dif- 
ferent products of spins residing on four neighboring sites, 
including, e.g., the combination (7„(7n+3- The number of 
quasi- shift operators will be also three instead of the 




FIG. 2: Chiral-spin order parameter, versus k — " ^ . 



two as in Eq. (|10p . It is plausible to think that one can 
use relations analogous to Eqs. (|18p - ((20|) in order to 
extract contributions of different summands out of the 
local Hamiltonian operators, at least in the homogeneous 
limit, LUi = UJ2 = 0. These investigations may constitute 
separate article. 
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